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In spite of the massive interest that the generalized Manakov equation has attracted in the past two decades, 
no physical system which is quantitatively described by this equation has been reported so far. In this paper 
we show that propagation in a group of degenerate modes of a multi-mode optical fiber satisfies this equation 
in the presence of random mode coupling. Consequently, this is the first reported physical system that admits 
true multi-component soliton solutions. The reported formalism constitutes the starting point for future studies 
of nonlinear effects in multi-mode fiber transmission. 
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Several years after the publication of the famous Za- 
kharov and Shabat paper [1] which showed that the 
scalar nonlinear Schrodinger equation (NLSE) is inte- 
grable, Manakov [2] discovered that a special form of 
a two-component vector NLSE in which the nonlinear- 
ity is isotropic shares the same property. However, no 
practical physical system was described by this equation 
until Wai, Menyuk and Chen [3] found that Manakov's 
equation accurately describes the propagation of a po- 
larized optical field in single-mode birefringent optical 
fibers in the presence of random mode coupling. This 
implied that such fibers support the existence of vector 
solitons, as predicted by Manakov almost two decades 
earlier. 

Recently, with fiber-communications exhausting the 
capacity of single-mode fibers, multi-mode fibers are 
being considered for transmission, with the purpose of 
increasing the overall fiber-communications throughput 
through spatial multiplexing. This new paradigm calls 
for the extension of the existing theoretical framework to 
the multi-mode fiber case. Random coupling of degener- 
ate or quasi degenerate modes is a distinctive feature of 
these fibers [4, 5]. We have addressed the linear problem 
of mode coupling and modal dispersion in multi-mode 
fibers in a recent paper [6] . In this paper we focus on non- 
linear propagation regime, where we show that the evo- 
lution of the electric field in a group of degenerate modes 
of a multi-mode optical fiber is described by a general- 
ized multi-component Manakov equation [7] . The central 
feature of the generalized Manakov equation is that it is 
integrable [7] and hence supports the propagation of soli- 
tons. Those are particle-like waveforms that remain un- 
changed in the process of propagation and preserve their 
identity when colliding with each other. It is the latter 
property that distinguishes between solitons and generic 
solitary waves, whose shape is also unaltered by propa- 
gation, but becomes corrupted when two, or more, such 
waveforms collide. Interestingly, solitary waves in cer- 
tain multi-mode optical systems were shown to be pos- 
sible using specific combinations of physical parameters 
[8, 9, 10, 11, 12]. Yet, unlike the case which wc present 



here, these systems are modeled by non-integrable equa- 
tions and hence they do not support true soliton propa- 
gation. 

We consider propagation over a group of degenerate 
modes in a multi-mode optical fiber. The degenaracy of 
the modes implies that fiber imperfections in the form 
of mechanical stress or manufacturing distortions pro- 
duce strong random coupling between them on a length- 
scale typically much shorter than the length-scale of the 
nonlinear evolution [13]. On the other hand, the distinc- 
tive difference in wave vector between non degenerate 
groups of modes makes coupling between them signifi- 
cantly smaller and hence we neglect it in our analysis. 
The four degenerate LPn modes [4] of a step- index fiber 
are a relevant example for a group of degenerate modes 
of the kind that wc consider here. Larger groups of de- 
generate modes can be encountered in situations where 
multiple multi-mode fibers are combined in a multi-core 
fiber structure [14]. 

Starting from the coupled NLSE and assuming ran- 
dom mode coupling, we generalize the standard Manakov 
equation [3] to the multi-mode case. We verify the accu- 
racy of the generalized Manakov equation and demon- 
strate the existence of multidimansional vector solitons, 
by solving the complete set of coupled NLSE numerically. 

The electric field in a group of A^ degenerate spatial 
modes is represented by a 2iV-dimensional complex val- 
ued vector E{z,t), which is constructed by stacking the 
Jones vectors of the N individual spatial modes one on 
top of the other. The term spatial mode is used here and 
throughout the paper to refer to the set of two polariza- 
tion modes sharing, in the weakly guiding fiber approx- 
imation [15], the same lateral field profile. The compo- 
nents of E{z, t) then satisfy the following set of coupled 
NLSE [16, 17] 



dz 



^E 



,(1) 



Er, 



^jm-^m 



T.p'> 



(1) dE„ 



dt 



m hkm 



which, in a vector form, can be re-expressed as 
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where B^'^ {z),i = 0,1,2 are 2N x 2N Hermitian matrices 
and where we use Bj with j ~ 1, • ■ • , 'ZN to denote the 
set of complex orthogonal unit vectors used to represent 
the electric field. 

The fourth term on the right-hand side of Eqs. (1) 
and (2) represents the coupling induced by the Kerr non- 
linearity of the fiber. The coefficient 7 = Ldon2 / cA^s is 
identical to the usual nonlinearity coefficient appearing 
in the scalar NLSE [16] of a single mode fiber, where 712 
is the Kerr coefficient of glass, c is the speed of light in 
vacuum, and Acs is the effective area of the fundamental 
mode at central frequency ujq. The dimcnsionlcss con- 
stants Cjhkm depend on the details of the spatial mode 
profiles and are obtained as follows [17]. Defining 
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where Fn{x,y) describes the lateral mode profiles, and 
nf{x,y) is the refractive index, we obtain 
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[2D\\[,+D\',[,]/3 = D[{[„ so that din 
C2222 = 1, where the indices 1 and 2 are used to denote 
the two polarizations of the fundamental mode. Note 



that ^0 



cAcs, where rics is the effective refractive 



index of the fundamental mode at the central frequency 
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The first three terms on the right-hand side of Eqs. (1) 
and (2) account for linear propagation and coupling [6]. 
While in the ideal case the matrices B(*)(z) {i = 0, 1, 2) 
are proportional to the identity, unavoidable position 
dependent perturbations result in the presence of off- 
diagonal terms, producing linear coupling between the 
various modes. The strongest coupling results from the 
term proportional to 'B^^' , and in most cases of practi- 
cal interest, its characteristic length-scale is shorter by 
orders of magnitude than the length-scale that character- 
izes the nonlinear evolution [16]. Under these conditions 
the orientation of the electric field vector, defined by 
E/\E\, uniformly samples the space of possible orienta- 
tions within the length-scale of the nonlinear evolution. 
Hence, field propagation can be well approximated by 
averaging the nonlinear terms with respect to the elec- 
tric field's orientation. As can be anticipated based on 



symmetry and dimensionality arguments, the averaged 
nonlinear term of (2) must reduce to the form 
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where k is a dimensionless parameter that depends only 
on the nonlinear coupling coefficients Cjhkm- Following 
the steps described in the appendix, we show that 
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where Sij is Kronecker's delta function. Equation (2) can 
thus be reduced to 

This equation is a 2N component generalization of 
the well-known Manakov-PMD equation describing the 
propagation in a single-mode optical fiber with random 
polarization coupling [18]. Indeed, when A^ = 1, we have 
C1221 = C2112 = 2/3 and Cnn = C2222 = 1 [16], yield- 
ing K = 8/9, as expected [18]. The non-diagonal terms 
in the matrices B^^^^ and B^^^ are due to the frequency 
dependence of the mode coupling, which is typically only 
a small correction to the coupling contained in the fre- 
quency independent part B(°) . By neglecting these terms 
and transforming Eq. (9) to a reference frame evolving 
with B'^*'), the final form of the generalized Manakov 
equation is obtained 
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The terms /3' and /3" are the inverse group velocity and 
the dispersion coefficient, respectively. The form of Eq. 
(10) is identical to that of the familiar Manakov equation 
[2], except that it describes the propagation of signals 
through 2N degenerate propagation modes of a multi- 
mode fiber. 

The generalized Manakov equation (10) is integrable 
by the inverse scattering transform [7], and hence it ad- 
nfits the propagation of solitons. Unlike generic solitary 
waveforms that were shown to exist in multi-mode fibers 
with certain parameter combinations [8], true solitons 
can only exist in the strong coupling regime, where prop- 
agation is accurately described by Eq. (10). Thus, in or- 
der to test the accuracy of Eq. (10) in describing multi- 
mode propagation, we verify that the fundamental soli- 
ton waveform, which is a rigorous analytical solutions 
of Eq. (10), indeed forms a solitary solution when the 
coupled NLSE (1) are solved numerically in the strong 
mode-coupling regime. In addition, we check that when 
two such waveforms are launched into different fiber 
modes and at different central optical frequencies, they 
collide elastically without leaving a trace in the form of 
a dispersive wave. As noted earlier, this phenomenon is 
a distinctive feature of solitons [1] distinguishing them 



from generic solitary waves. We find that Manakov soli- 
tons can be observed when the length-scale characteriz- 
ing the correlation of mode coupling is smaller than the 
soliton length Lg = t^/|/3"| (related to the soliton period 
Zs by Zs = iiLs/2). 

In the simulations we consider propagation in the two 
degenerate LPn modes of a step-index optical fiber, 
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LP^'J'' and IjP\'{ [19]. We used a core radius of 7.5/im, 
a core refractive index of 1.4621, a refractive index 
step of 9.7 X 10"'^ and a dispersion coefficient /?" = 
— 25ps^/kni. The nonlinear coefficient in our computa- 
tions was 712 = 2.6 X 10^^°m^W~^, corresponding to 



^km- 



0.76, and the effective area 



7 ~ 0.835 W- 

was of 126/ini^ for the fundamental mode. 

In Fig. 1 we consider the case in which we launch 
the waveform describing the fundamental soliton of 

1 /2 

the Manakov equation (10), Pq sech(t/T) with Pq = 
/3"|/(7KT^), into each of the two modes LP^°^ and LP^^\ 
The two launched soliton waveforms are separated by 
T — 30r in time and by Aw = 0.4/r in frequency. As ran- 
dom coupling of polarizations within each spatial mode is 
always very fast [3] , the launch polarization states of the 
two waveforms is immaterial. We assume a soliton half- 
width Th ~ 17ps corresponding to r = 30ps. The soliton 
length is Lg — 36km and the peak power Po = 44mW. 
The correlation length of the mode coupling was taken 
to be ic = 100m [Lc/Lg ~ 2.8 x 10"^). It is evident 
that the launched waveforms are indeed solitary solu- 
tions in this regime and the collision between them is 
elastic, as expected from true soliton behavior. The col- 
lision dynamics shown in Fig. 1 is indistinguishable from 
that obtained integrating the Manakov equation with the 
same initial conditions. 

In Fig. 2 we consider the regime in which the spatial 
modes are perfectly isolated so that no coupling between 
them occurs, although full coupling occurs between po- 
larizations in each mode. In this case we launch the same 
scch-waveform, but with Pq = |/3"|/(87C33337r^/9), 
where C3333 ~ 1.07 corresponds to self- phase modula- 
tion in the LPn modes (see definitions around Eq. (6), 
where 3 is the index of mode LP]^" ). Since the launched 
waveform is solitary for the individual modes LPj^° and 
LP^]^^, the pulses propagate unperturbed until their colli- 
sion. Yet, upon collision they stick to each other, forming 
a new (and non-solitary) pulse, while some energy is lost 
in the form of dispersive radiation. This behavior is in 
clear contrast with the elastic collision observed in Fig. 
1, which corresponds to the strong coupling regime, in 
which the generalized Manakov equation holds. 

In Figs. 3 and 4 we consider cases where the correla- 
tion length of the random mode coupling is L^ ~ 10km 
{L^/Ls ~ 2.8x10-1) and Lc = 100km (Lc/L^ ~ 2.8), re- 
spectively. In the case of L^ = 10km shown in Fig. 3, the 
pulse evolution is similar to that observed in the strong 
coupling case of Fig. 1, although some perturbations to 
the propagating waveforms can be observed due to the 
insufficient averaging of the random mode coupling in 
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Fig. 1. Collision of two solitons initially in two orthog- 
onal modes in the case of full coupling, and initial an- 
gular frequency separation Ao; = 0.4/t. The time axis 
is normalized to the soliton time r, the spatial axis to 
the soliton length Lg and the vertical axis to the soliton 
peak power P^. The ratio of the correlation length to the 
soliton length is Lc/Lg — 2.8 x lO"'^. 
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Fig. 2. Collision of two solitons initially in two orthogo- 
nal modes in the case of absence of coupling, and initial 
angular frequency separation Aw = 0.4/t. The time axis 
is normalized to the soliton time r, the spatial axis to 
the soliton length Lg and the vertical axis to the soliton 
peak power P„ 



this regime. These perturbations are further exacerbated 
when the correlation length of the mode coupling is ex- 
tended to 100km, as shown in Fig. 4. In both Figs. 3 
and 4, the launched waveforms were Manakov-equation 
solitons, identical to those used in Fig. 1. 

The generalized Manakov equation (10) constitutes a 
powerful tool for understanding the nonlinear dynam- 
ics in multi-mode optical fibers. As in the case of its 
two-component predecessor, describing polarization dy- 
namics in single- mode fibers [3], it constitutes a start- 
ing point for most analytical descriptions of nonlinear 
propagation phenomena [20] . In addition, as it incorpo- 
rates the effect of random coupling analytically, it can 
be used for an effective simplification of numerical stud- 
ies of signal propagation [21]. Possible applications in- 
clude the analysis of inter and intra-channel interfer- 



\E(z,t)\ /P„ 




Fig. 3. Collision of two solitons initially in two orthogonal 
modes in the same condition of Fig. 1, with ratio of the 
correlation length to the soliton length of Lc/Lg — 0.28 x 
10-1. 
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Fig. 4. Collision of two solitons initially in two orthogonal 
modes in the same condition of Fig. 1, with ratio of the 
correlation length to the soliton length of Lc/Ls — 2.8. 



ence in fiber-communications systems using multiplexing 
in both wavelength and space. Recently, such systems, 
deployed over multi-mode [22] and multi-core [23] op- 
tical fibers are being hectically studied with the hope 
of providing an economically viable path for overcom- 
ing the capacity crunch experienced in the field of fiber- 
communications [24]. 

In order to show that in the presence of random mode 
coupling the coupled NLSE (1) are equivalent to the 
generalized Manakov equation (10) we solved (1) nu- 
merically and showed that it supported the fundamental 
soliton waveforms, which are exact analytical solutions 
of (10). To demonstrate that those were indeed solitons 
and not generic solitary waves, we verified that collisions 
between pairs of such waveforms were elastic, namely 
such that maintain the identity of the individual soli- 
tons and produce no dispersive waves. The accuracy of 
the Manakov equation in describing the field evolution 
was shown to be excellent when the correlation length of 
the mode coupling was smaller than the nonlinear length 
by an order of magnitude, or more. 

The fact that light propagating through degenerate 
modes in a multi-mode fiber satisfies the generalized 



Manakov equation can be exploited for the prediction 
of multiple general properties of signal propagation. An 
important example is that the isotropy of the Manakov 
equation implies that collisions between two polarized 
pulses (i.e. pulses whose orientation E{t)/\E{t)\ is time 
independent) can always be considered as occurring in a 
two dimensional (complex) subspace of the 2N dimen- 
sional space spanned by the entire set of modes [25, 26]. 
Consequently, the description of two-pulse interactions 
becomes identical to the description of two-pulse inter- 
actions in single-mode fibers, allowing use of the conven- 
tional Stokes-space representation [27]. Then, it can be 
shown that pulse collisions correspond to the Stokes vec- 
tor's representing each of the pulse precessing about each 
other [20, 28]. When the colliding pulses are orthogonal 
to each other (implying that their corresponding Stokes 
vectors are antiparallel to one another) their orientations 
in the 2A'^-dimensional space of the electric field remain 
unaltered. Conversely, when the input pulses are not or- 
thogonal, their Stokes vectors are not aligned with each 
other and precession in Stokes space implies that the 
pulses' orientations are modified by the collision, so that 
the distribution of the overall optical power among the 
various fiber modes changes. This effect, which could not 
be easily predicted using Eq. (1), may be of importance 
to space and wavelength multiplexed systems when op- 
tical nonlinearities are non negligible. 

We showed that propagation in a degenerate group of 
randomly coupled spatial modes of a multi-mode opti- 
cal fiber is the first reported physical scenario which is 
described by the generalization [7] of the famous Man- 
akov equation [2], and hence admits true vector soliton 
solutions. The key parameter of the equation describing 
the nonlinear effect was expressed in terms of the stan- 
dard parameters of a generic optical fiber. This equation 
constitutes a starting point for analytical and numerical 
studies of nonlinear effects in multi-mode transmission. 
Its importance is emphasized by the massive recent in- 
terest in spatially multiplexed transmission using multi- 
mode and multi-core optical fibers. 

Appendix 

In order to derive Eq. (8), we perform scalar multipli- 
cation of both sides of Eq. (7) by E and average with 
respect to the field's orientation. This yields 



jhkm 



(11) 



where E denotes statistical averaging and Qjhkm = 
E [E'^EkEniE*] /\E\'^. Since E can be considered as a 
constant modulus vector whose orientation is uniformly 
distributed, the statistical average can be performed as 
follows. We introduce an auxiliary random vector X 
with 2A^ complex, statistically independent components. 
The real and imaginary parts of each component Xi of 
X are statistically independent standard Gaussian vari- 
ables having zero-mean and unit variance. It can be ar- 
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gued that Qjhk^n = E 

tiplying both sides by X|* and performing another av- 
erage (with respect to the square modulus |^P), we find 
that 
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Since ah components of X are statistically independent 
standard complex Gaussian variables, the numerator is 
4 {Shk^jm + ^hmSjk)- The denominator is the second mo- 
ment of a chi-square random variable having AN degrees 
of freedom and hence its value is 4A^(4A^-|-2). Using these 
results, the final form of Eq. (8) readily follows. 
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